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1.  SUMMARY 

In  the  second  half  of  the  third  year  of  the  project,  further  experiments  have  been 
conducted  in  several  areas.  These  include  detailed  PIV  and  LDV  measurements  of  the 
vortical  flow  in  the  wind  tunnel  to  provide  high  Reynolds  number  data  for  comparisons 
with  the  computational  simulations  (carried  out  by  Dr.  Ray  Gordnier  and  Dr.  Miguel 
Visbal,  Computational  Sciences  Branch,  Air  Vehicles  Directorate),  hot-wire 
measurements  for  spectral  analysis  of  unsteady  flow  features,  as  well  as  flow 
visualization  for  oscillating  rigid  wings  in  the  water  tunnel  to  study  the  effects  of  wing 
sweep  angle  and  symmetric  perturbations.  In  both  the  wind  tunnel  and  water  tunnel 
experiments,  flow  visualization,  velocity  measurements  and  force  measurements  provide 
insight  to  lift  enhancement  and  flow  reattachment  over  flexible  and  oscillating  wings.  A 
list  of  publications  originated  from  this  project  and  two  recent  conference  papers  are 
attached. 


2.  MAIN  FINDINGS 

2.1.  Wind  tunnel  experiments 

2.1.1.  PIV  and  LDV  measurements 

The  main  objective  of  these  measurements  was  to  provide  high  Reynolds  number 
data  for  the  vortical  flow  before  the  stall  (a=15°).  These  data  have  been  used  for 
comparisons  with  the  computational  simulations.  Figure  1  shows  the  PIV  data  from  the 
wind  tunnel  experiments  and  also  compares  with  the  water  tunnel  experiments.  There  are 
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strong  similarities,  although  some  differences  in  the  breakdown  location  and  the  spanwise 
location  of  the  vortex  exist.  LDV  measurements  provide  higher  resolution  and  also 
provide  information  on  the  rms  velocity.  An  example  is  shown  in  Figure  2.  Detailed 
results  (not  shown  here)  were  provided  to  our  collaborators  in  the  AFRL. 

2,1.2.  Spectral  analysis 

Both  LDV  and  hot-wire  measurements  were  conducted  for  spectral  analysis  of 
unsteady  flow  features.  The  data  rate  of  the  LDV  data  was  not  sufficiently  high  for  this 
purpose.  Therefore,  extensive  measurements  were  taken  with  a  single  hot  wire,  and  also 
with  two  hot-wires  (for  simultaneous  measurements).  Figure  3  shows  the  spectra  of 
velocity  fluctuations  at  different  spanwise  locations  for  the  rigid  wing,  a=20°,  x/c=0.7. 
For  this  angle  of  attack,  vortex  breakdown  is  already  at  the  apex  of  the  wing.  It  is  seen 
that  a  broad  peak  exists  in  the  spectra,  and  also  the  center  frequency  is  approximately  the 
same.  As  the  locations  of  measurements  are  close  to  the  separated  shear  layer,  it  is 
believed  that  these  peaks  correspond  to  the  dominant  frequency  of  the  shear  layer 
structures. 

Figure  4  shows  the  effect  of  streamwise  distance  on  the  spectra  of  velocity 
flucutuations  for  a=20°.  It  is  seen  that  the  center  peaks  decrease  with  increasing 
streamwise  distance.  Figure  5  shows  the  spectra  for  various  angles  of  attack  in  the  range 
of  a=5°  to  a=25°  for  the  rigid  wing.  The  dominant  peaks  are  seen  for  a  >  15°,  which 
covers  the  flow  regimes  when  vortex  breakdown  exist  over  the  wing  as  well  as  stalled 
flow.  Hence,  it  is  believed  that  these  are  due  to  the  shear  layer  instabilities,  although 
there  is  the  possibility  that  vortex  breakdown  might  be  driving  these  instabilities  at  lower 
incidences. 

Simultaneous  measurements  at  two  locations,  which  are  the  mirror  images, 
indicated  that  there  was  no  coherence  between  the  two  signals  as  shown  in  Figure  6.  This 
was  somewhat  surprising  in  that  the  dominant  peaks  observed  in  the  spectra  are  not 
correlated  at  all  for  the  rigid  wing.  When  similar  measurements  were  carried  out  for  the 
flexible  wing,  the  large  sharp  peak  corresponding  to  the  wing  vibrations  is  seen  in  the 
spectrum  of  velocity  fluctuations  in  Figure  7  in  the  post-stall  region.  Figure  8  shows  the 
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amplitude  and  phase  of  the  coherence  between  the  two  velocity  signals,  which  confirm 
the  anti- symmetric  nature  of  the  oscillations. 

2.2.  Water  tunnel  experiments 

2.2.1.  Laser  fluorescence  flow  visualization  for  oscillating  wings 

Experiments  with  a  rigid  wing  undergoing  small-amplitude  rolling  motion  were 
conducted  in  order  to  simulate  the  effect  of  antisymmetric  vibrations.  In  these 
experiments,  fluorescent  dye  was  released  from  a  slot  near  the  leading-edge,  hence  it 
“marks”  the  vorticity  shed  with  separation.  The  objective  was  to  visualize  the  shear  layer 
structures. 

Figure  9  shows  flow  visualization  for  the  stationary  and  rolling  delta  wing  (with 
an  amplitude  of  1°  and  St=1.0)  for  angles  of  attack  of  a=15°,  a=20°  and  a=25°  at 
x/c=0.8.  While  the  difference  is  small  for  a=15°,  the  change  in  the  size  of  the  separation 
region  is  considerable  for  the  oscillating  wing  due  to  earlier  reattachment.  The 
wavelength  of  the  vortical  structures  on  stationary  and  oscillating  wing  suggests  that  the 
natural  frequency  of  the  shear  layer  instability  is  close  to  St=l,  which  is  also  close  to  the 
center  peak  frequency  of  the  hot-wire  measurements.  Recall  that  this  also  corresponds  to 
an  “optimum  frequency”  for  which  breakdown  location  is  delayed  to  a  maximum 
distance  from  the  apex.  Figure  10  shows  the  time-averaged  flow  visualization  for  the 
three  incidences. 

2.2.2.  Effect  of  wing  sweep 

Figure  1 1  and  12  show  flow  visualization  pictures  for  A=40°  and  A=30°  wings  at 
various  post-stall  angles  of  attack  for  stationary  and  oscillating  cases.  It  is  seen  that  in  all 
cases  flow  reattachment  occurs  over  the  wing  for  the  roll  oscillations  at  St=1.0.  Although 
not  shown  here,  detailed  examination  and  measurement  of  reattachment  region  showed 
that  St=l  is  around  the  “optimum”  frequency,  which  is  similar  to  A=50°  wing  results. 

2.2.3.  Symmetric  oscillations 

Effect  of  symmetric  perturbations  in  the  form  of  small  amplitude  pitching 
oscillations  (1°  amplitude)  were  studied.  The  results  of  flow  visualization  are  shown  in 
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Figure  13.  Again,  earlier  reattachment  and  even  the  re-fonnation  of  the  leading-edge 
vortex  is  observed  for  the  oscillating  wing.  Hence,  for  active  control  purposes,  both 
symmetric  and  anti-symmetric  forcing  would  work.  However,  passive  control  for  a 
flexible  wing  occurs  only  in  the  anti-symmetric  mode. 
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(b)  a  =  15°,  <)>  =  10°,  Rec  =  620,000,  Wind-tunnel. 


Figure  1:  Comparison  of  PIV  measurements  for  wind  and  water  tunnel  experiments. 
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Horizontal  traverse  at  x/c=20%,  a=15°,  z/s=0.17 


Figure  2:  LDV  measurements  of  velocity  profile  upstream  of  breakdown  for  a=15°. 
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Stationary  wing 


Oscillating  wing 


a=15° 


a=20° 


a=25° 


Figure  9:  Laser  fluorescence  flow  visualization  for  stationary  and  rolling  wings  (St=  1 .0,  A(j)=l°). 
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Stationary  wing 


Oscillating  wing 


Figure  10:  Time-averaged  laser  fluorescence  flow  visualization  for  stationary  and  rolling 
wings  (St=1.0,  A(j)=l°). 
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a  =  20° 


a  =  25° 


Figure  11:  Wing  sweep  angle,  A=40°:  flow  visualization  for  stationary  and  oscillating  wings 
in  roll  (St=1.0,  A(j)=5°). 
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a  =  15° 


a  =  20° 


Figure  12:  Wing  sweep  angle,  A=30°:  flow  visualization  for  stationary  and  oscillating  wings 
in  roll  (St=1.0,  A(j)=5°). 
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St  =  0 


St  =  0.5 


St  =  1.19 


Figure  13:  Effect  of  symmetric  (pitching)  oscillations:  flow  visualization  for  stationary  and 
oscillating  wings  (St=1.0,  Acj)=l0,  wing  sweep  angle  A=50°). 
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Passive  Flow  Control  over  Flexible  Non-Slender  Delta  Wings 

Gordon  Taylor*,  Andreas  Kroker1  and  Ismet  Gursuf 
University  of  Bath,  Bath,  BA2  7 AY,  United  Kingdom. 


Force  measurements  over  a  range  of  non-slender  delta  wings  have  demonstrated  the 
ability  of  a  flexible  wing  to  enhance  lift  and  delay  stall  compared  with  a  rigid  wing  of  similar 
geometry.  The  work  has  extended  the  results  of  a  recent  study  to  include  delta  wings  with  a 
range  of  sweep  angles.  It  has  been  shown  that  the  greatest  lift  enhancement  was  observed 
over  the  wing  with  smallest  leading  edge  sweep.  Additionally  for  the  only  slender  wing 
considered  no  lift  enhancement  was  observed  indicating  that  it  is  a  feature  of  the 
fundamentally  different  flow  that  occurs  over  low-sweep  wings  that  is  responsible  for  the 
phenomenon.  The  variation  of  RMS  lift  force  and  rolling  moment  coefficients  with 
incidence  for  all  the  wings  concerned  suggest  that  the  wings  are  vibrating  in  an  anti¬ 
symmetric  structural  mode  in  the  lift  enhancement  region.  This  supports  previous  evidence 
for  the  50°  wing  that  suggests  the  same.  Further,  experiments  of  a  half-wing  model  suggest 
that  this  wing  does  not  undergo  lift  enhancement,  and  as  such  it  may  be  that  the  anti¬ 
symmetric  vibration  is  a  necessary  condition  for  the  lift  enhancement  to  exist.  Wing-tip 
deflections  for  some  of  the  wings  were  also  measured,  which  showed  that  the  slender  wing 
undergoes  much  higher  time-averaged  deflections  than  the  non-slender  wings,  although  it 
also  experiences  a  much  lower  level  of  buffet.  PIV  and  LDV  measurements  have 
demonstrated  the  striking  difference  between  the  surface  flows  over  the  flexible  and  rigid 
wings  in  the  post-stall  region  at  the  same  incidence.  Implementing  flexibility  on  a  low  sweep 
wing  results  in  continued  reattachment  of  the  shear  layer  to  much  higher  incidences  than 
would  otherwise  be  observed  over  a  rigid  wing. 
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Nomenclature 

=  root  chord  length 
=  frequency 
=  local  semi-span 
=  wing  thickness 
=  lift  coefficient 
=  Rolling  moment  coefficient 
=  elastic  modulus 
=  chord  Reynolds  number 
=  Strouhal  number  (  =  fc/Uf) 

=  free-stream  velocity 
=  angle  of  attack 

=  wing-tip  vertical  deflection  under  load 

=  span-wise  bending  stiffness  normalized  by  free-stream  dynamic  pressure  (  =  12(1- V)pjjjs3/Et3) 
=  Poisson’s  ratio 
=  density 
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I.  Introduction 


he  recent  interest  in  Micro- Air- Vehicles  (MAVs)  and  Unmanned-Air- Vehicles  (UAVs)  has  resulted  in  a  need  to 

T further  our  understanding  of  flows  over  non-slender  delta  wing  configurations.  While  a  non-slender  wing  may  be 
out-performed  by  a  slender  wing  in  terms  of  maximum  achievable  lift  and  stall  angle1’2’3,  the  non-slender  wing 
offers  alternative  advantages  such  as  lower  lift-to-drag  ratio4  and  a  greater  surface  area  allowing  greater  options 
for  the  internal  packaging  of  flight  systems,  thus  minimizing  the  impact  on  external  aerodynamics.  Both  of  these 
benefits  are  particularly  important  for  small  aircraft  such  as  MAVs.  For  the  purposes  of  this  investigation,  a  non- 
slender  wing  is  defined  as  one  with  leading  edge  sweep  less  than  60°.  The  flows  that  occur  over  non-slender  wings 
have  been  shown  to  differ  substantially  from  those  documented  over  more  slender  planforms.  While  a  slender  wing 
produces  a  single  coherent  primary  vortex  from  each  leading  edge,  recent  computational5  and  experimental6  studies 
have  demonstrated  the  existence  of  a  ‘dual’  primary  vortex  structure  over  a  non-slender  wing.  This  vortex  structure 
is  a  direct  result  of  the  proximity  of  the  vortex  formation  to  the  wing  surface7'8,  which  results  in  significant 
interactions  between  the  secondary  and  primary  vortices.  The  primary  vortex  is  split  into  two  separate  vortices  by 
the  secondary  vortex,  which  occurs  much  further  away  from  the  primary  vortex  in  slender  wing  flows.  A  further 
consequence  of  the  interaction  between  the  boundary  layer  and  vortex  flows  is  the  sensitivity  of  non-slender  wing 
flows  to  Reynolds  number6  that  is  not  observed  over  slender  wings9.  For  Re  <  30,000,  decreasing  Reynolds 
number  results  in  a  downstream  movement  of  the  vortex  breakdown  and  a  displacement  of  the  vortex  core  towards 
the  wing  center-line6.  The  vortices  generated  over  low  sweep  delta  wing  undergo  breakdown  at  very  low 
incidences6;  for  a  50°  wing  the  breakdown  has  been  observed  over  the  upstream  half  of  the  wing  at  an  incidence  of 
just  2.5°.  Detailed  three-component  PIV  measurements  have  demonstrated  the  wake-like  velocity  profile  of  the 
leading  edge  vortices  produced  over  a  50°  wing10  at  low  Reynolds  numbers.  Further  measurements  have  shown  that 
this  is  a  feature  observed  at  low  incidences11,  and  is  possibly  a  result  of  the  propinquity  of  the  vortices  to  the  wing 
surface. 


While  the  flows  over  non-slender  wings  have  become  a  more  popular  topic  in  the  literature  over  recent  years, 
there  have  been  few  investigations  of  the  effects  of  flexibility.  An  initial  study  into  flexibility  of  a  60°  wing  showed 
that  significant  buffeting  was  present  in  the  region  just  prior  to  the  stall12.  For  this  wing,  of  thickness-to-chord  ratio, 
t/c  =  0.44,  at  a  Reynolds  number,  Re  =  O(106),  finite  element  analysis  complemented  with  wing-tip  accelerometer 
measurements  indicated  that  in  this  state  the  wing  was  vibrating  in  the  second  anti-symmetric  structural  mode12. 
This  was  later  confirmed  with  computational  simulations,  which  showed  significant  peaks  in  the  spectra  of  pressure 
fluctuations  centered  on  the  2nd  and  3rd  modes  of  vibration13.  Only  recently  have  the  effects  of  flexibility  been 
studied  for  a  non-slender  delta  wing.  Studies  of  a  50°  wing  of  thickness,  t/c  =  0.65%  at  a  similar  Reynolds  number 
have  also  shown  significant  levels  of  buffet11.  There  are  a  number  of  sources  of  unsteadiness  that  may  manifest 
themselves  as  buffet,  and  these  sources  have  been  studied  extensively  over  slender  delta  wings14’15.  It  has  been 
shown  that  for  a  slender  wing  the  onset  of  breakdown  forward  of  the  trailing  edge  results  in  a  sharp  increase  in  the 
level  of  buffet,  as  measured  by  normal  force  fluctuations1  and  RMS  wing-root  strain16.  Further,  the  onset  of  stall 
corresponds  to  a  drop  in  the  level  of  buffet12,  implying  that  the  instabilities  associated  with  the  breakdown  process 
are  the  dominating  influence  in  the  buffeting  process  over  slender  wings17.  Flowever,  over  non-slender  wings, 
breakdown  of  the  vortices  has  been  observed  over  the  wing  at  incidences  as  low  as  2.5°,  and  breakdown  can  reach 
the  apex  at  an  incidence  much  lower  than  the  stall  angle6.  In  spite  of  this  behavior,  peak  buffeting  still  occurs 
around  the  stall  the  angle  over  low-sweep  wings11,  and  instabilities  associated  with  the  shear  layer  reattachment 
process  at  high  incidences  have  been  proposed  as  the  leading  cause  of  buffet  over  these  wings. 

A  recent  study  has  documented  a  particularly  interesting  feature  of  non-slender  wings  of  sufficient  flexibility15. 
It  was  shown  that  for  a  wing  of  thickness,  t/c  =  0.32%  at  a  Reynolds  number,  Re  =  6.2e5,  a  region  of  significant  lift 
enhancement  exists  in  the  range  of  incidences  immediately  following  the  stall.  In  this  region,  lift  was  shown  to 
increase  by  45%  at  a  =  30°,  and  the  stall  may  be  delayed  by  up  to  9°.  Large  wing  tip  deflections,  both  mean  and 
fluctuating,  were  associated  with  the  lift  enhancement  region,  along  with  a  switching  of  the  dominant  structural 
mode  from  the  fundamental  to  the  second  anti-symmetric.  The  large  time-averaged  deflection  was  initially 
suspected  as  the  cause,  but  comparison  with  wings  with  span- wise  camber19  indicated  that  this  should,  if  anything, 
result  in  a  reduction  in  lift.  Detailed  measurements18  indicate  that  it  was  the  magnitude  of  the  tip  vibrations  that 
were  responsible  for  this  phenomenon. 

The  aim  of  this  investigation  was  to  study  further  the  lift  enhancement  phenomenon  that  has  been  observed  over 
non-slender  delta  wings.  Specifically,  the  aims  were  to:  (1)  determine  whether  the  same  lift  enhancement  could  be 
achieved  over  wings  with  leading  edge  sweep  angles  other  than  50°;  and  (2)  undertake  measurements  to  document 
the  flow  structure  in  the  lift  enhancement  region.  These  aims  were  achieved  by  complementing  lift  force  data  with 
near-surface  PIV  and  off-surface  LDV  measurements. 
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II.  Methodology 

Wind  tunnel  experiments  were  conducted  in  the  high  speed  working  section  of  the  7'  x  5'  closed  circuit  facility  at 
the  University  of  Bath.  Models  were  supported  using  the  high-incidence  mechanism  as  shown  in  Figure  1. 
Experiments  were  conducted  at  a  free  stream  velocity  of  Uoo  =  31.7  m/s.  Maximum  blockage  for  the  wind  tunnel 
was  approximately  2.3%  for  the  60°  wing  at  a  =  50°.  Models  of  varying  leading  edge  sweep  and  thickness  were 
tested.  Each  model  was  a  planar  delta  wing  with  a  pressure  surface  leading  edge  bevel  of  45°,  and  a  square  trailing 
edge.  Figure  3  shows  a  schematic  of  the  wing  of  50°  leading  edge  sweep;  in  addition  to  this,  wings  of  40°,  45°,  55°, 
and  60°  sweep  were  tested.  Each  model  had  a  semi-span,  s  =  0.260  m  such  that  the  span-wise  flexural 
characteristics  of  the  thin  wings  remained  consistent  between  the  models.  Chord  lengths,  c,  in  the  range  0.218  - 
0.450  m  were  used  resulting  in  chord  Reynolds  numbers  in  the  range  440,000  -  900,000.  Each  model  was  made 
from  aluminium  alloy  1050A,  which  has  the  following  material  properties:  p  =  2700  Kg/m3;  E  =  69  GPa;  and  v  = 
0.3.  Models  of  thickness  1  mm  (t/s  =  0.38%)  and  5  mm  (t/s  =  1.92%)  were  tested.  The  1  mm  wing  was  designed  for 
high  flexibility,  and  achieved  significant  out-of-plane  deflections  as  documented  in  reference  18,  while  the  5  mm 
wing  was  rigid.  The  parameter  Xs,  the  reciprocal  of  the  spanwise  bending  stiffness  of  the  wing  normalized  by  the 
free-stream  dynamic  pressure,  was  used  to  give  an  indication  of  the  relative  flexibility  of  the  wings,  whereby  a 
higher  value  of  Xs  corresponds  to  a  more  flexible  wing.  The  above  conditions  resulted  in  a  non-dimensional  span- 
wise  bending  stiffness,  =  3.1  for  the  1  mm  wing  of  all  sweep  angles,  compared  with  a  value  of  =  0.025  for  the 
rigid  wings. 

Force  balance  measurements  were  undertaken  using  a  six-component  strain-gauged  internal  balance.  Forces 
were  normalised  by  qS',  moments  by  qSc,  where  q  is  the  free-stream  dynamic  pressure,  S  is  the  wing  reference  area, 
and  c  is  the  root  chord.  Root-Mean-Square  (RMS)  quantities  are  presented  following  filtering  of  the  voltage  signal 
to  remove  unwanted  electrical  noise.  Wing-tip  displacements  were  measured  using  a  high-speed  Kodak  digital 
camera  with  a  resolution  of  512  x  240  pixels,  and  a  capture  rate  of  500  frames  per  second.  Tip  displacements  were 
measured  from  the  resulting  images  and  calibrated  to  give  measurements  in  SI  units.  Buffet  response  was  measured 
by  mounting  Entran  EGA  miniature  accelerometers  with  a  mass  of  approximately  0.5  g  on  each  tip  of  the  flexible 
wings.  Signals  from  the  accelerometers  were  recorded  using  a  desktop  PC  with  a  12-bit  A/D  data  acquisition  card. 
Signals  were  recorded  for  a  period  of  5  s  at  a  sampling  frequency  of  1000  Flz,  resulting  in  a  record  spanning  over 
250  c/Uoo. 

Quantitative  flow  measurements  were  undertaken  using  PIV  and  LDV.  For  PIV,  illumination  of  the  plane  of 
interest  was  achieved  using  a  pair  of  pulsed  mini  Nd:YAG  lasers  with  a  maximum  energy  of  120  mJ  per  pulse. 
Images  were  captured  using  an  8-bit  digital  camera  with  a  resolution  of  4.2  million  pixels.  Each  case  was  averaged 
over  100  frames  at  a  frame  rate  of  3.75  Hz  in  order  to  yield  a  time-average  spanning  over  2500  c/U-,.,  and  velocity 
vectors  were  calculated  at  over  50,000  points  in  each  field.  The  flow  was  seeded  using  a  commercially  available 
theatre-smoke  generator.  The  measurement  planes  for  the  rigid  and  flexible  wings  are  shown  in  Figure  2(a).  In 
practice,  the  distance  from  the  wing  surface  was  prescribed  by  the  need  to  reduce  surface  reflections,  and  was  in 
most  cases  around  1mm  (h/c  =  0.003).  For  the  flexible  wing,  the  deformation  of  the  wing-tips  interfered  with  this 
measurement  plane,  and  so  a  plane  offset  from  that  on  which  the  apex  and  deformed  wing-tips  lay  was  chosen. 

LDV  measurements  were  undertaken  using  a  300  mW  air-cooled  Argon-Ion  laser  and  TSI  burst  correlator  unit. 
The  data  rate  was  of  the  order  of  200  Hz  with  a  burst  efficiency  of  approximately  50%,  and  a  total  of  5,000  data 
points  were  recorded  for  each  test  giving  a  sample  period  of  well  over  100  c/UM.  Measurements  of  the  axial  velocity 
were  taken  in  the  cross  flow  plane  at  x/c  =  0.4  and  0.7.  Single  traverses  vertically  above  the  centerline  were  also 
taken  at  a  higher  resolution,  and  the  resulting  velocity  profiles  for  the  rigid  and  flexible  wings  were  compared.  The 
measurement  planes  and  traverses  are  shown  in  Figure  2(b). 


III.  Results 


A.  Forces  and  Wing  Deformation 

Figure  4  compares  the  variation  of  lift  coefficient  with  incidence  for  the  flexible  and  rigid  wings  of  all  sweep 
angles  tested.  The  results  show  that  the  lift  enhancement  phenomenon  is  not  limited  to  wings  of  50°  leading  edge 
sweep;  rather,  the  effect  of  flexibility  was  to  result  in  a  region  of  enhanced  lift  for  all  wings  except  the  60°  wing. 
The  magnitude  and  extent  of  the  lift  enhancement  region  was  a  function  of  sweep  angle,  with  the  greatest 
enhancement  being  observed  over  the  wing  of  lowest  sweep.  This  is  further  demonstrated  in  Table  1,  which  lists 
some  statistics  regarding  the  lift  curves  shown  in  Figure  4.  In  this  table,  ACLmax  represents  the  increment  in 
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maximum  lift  coefficient  achieved  by  the  flexible  wing  compared  with  (a)  CLmax  of  the  rigid  wing,  and  (b)  Cl  of  the 
rigid  wing  at  (Xstall  of  the  flexible  wing. 

The  40°  wing  improved  the  maximum  lift  coefficient  by  over  40%,  and  at  the  point  of  stall  of  the  flexible  wing, 
lift  coefficient  was  enhanced  by  over  50%.  As  sweep  angle  was  increased  the  magnitude  of  the  lift  enhancement 
reduced,  with  the  overall  maximum  lift  coefficient  of  the  55°  wing  increased  by  just  1%,  while  the  60°  wing 
experienced  an  overall  reduction  in  lift.  Similarly,  for  the  40°  wing  a  delay  in  the  onset  of  stall  of  7°  was  achieved 
by  the  flexible  wing,  while  no  delay  in  stall  was  achieved  for  the  60°  wing.  For  all  the  wings,  the  nose-down 
pitching  moment  measured  about  the  apex  underwent  a  similar  magnitude  increase  in  the  lift  enhancement  region. 


C-Lmax 

ACLmax 

A 

Rigid 

Flexible 

(a) 

(b) 

Accstall 

40 

0.77 

1.09 

42% 

53% 

7° 

45 

0.83 

1.14 

38% 

49% 

6° 

50 

0.92 

1.13 

23% 

45% 

5° 

55 

1.07 

1.08 

1% 

23% 

3° 

60 

1.24 

1.18 

-5% 

-5% 

0 

Table  1:  Statistics  of  lift  enhancement  and  stall  delay  over  flexible  delta  wings  of  varying  sweep  angle.  ACLmax 
(a)  and  (b)  are  defined  in  the  text. 

Due  to  the  degree  of  flexibility  employed  in  these  experiments,  considerable  time-averaged  and  fluctuating 
displacements  were  observed.  For  the  wings  of  40°,  50°  and  60°  sweep,  the  variation  of  wing-tip  displacement  with 
incidence  is  shown  in  Figure  5.  The  data  for  the  50°  wing  was  previously  presented  in  reference  18.  For  the  40° 
and  50°  wings  the  lift  enhancement  region  was  accompanied  by  an  increase  in  both  the  amplitude  and,  to  a  lesser 
extent,  mean  (time-averaged)  tip  deflection.  For  the  60°  wing,  which  did  not  experience  an  increase  in  lift,  no  such 
discontinuity  in  the  displacement  curves  was  observed.  In  fact,  it  is  clear  that  the  60°  wing  experiences  much  less 
buffet  that  the  lower  sweep  wings,  although  the  peak  mean  deflection  is  greater.  Although  the  wings  were  carefully 
designed  to  have  identical  values  of  7,s,  and  therefore  equivalent  bending  stiffness,  the  chord  lengths  of  the  wings  are 
not  equal  resulting  in  a  lower  thickness-to-chord  ratio  for  the  60°  wing,  and  this  probably  accounts  for  the  greater 
span-wise  bending.  The  effect  of  a  span-wise  dihedral  camber  of  this  nature  is  to  reduce  the  effective  leading  edge 
sweep  angle  and  incidence  of  the  wing19,  and  these  effects  are  consistent  with  a  reduction  in  lift.  This  explains  why, 
in  the  absence  of  a  lift  enhancement  region,  the  flexible  60°  wing  underwent  a  reduction  in  the  lift-curve  slope 
compared  with  the  rigid  wing. 

Figure  6  shows  the  variation  of  RMS  lift  coefficient  with  incidence  for  all  the  wings.  While  there  is  clearly  an 
increase  in  the  level  of  buffet  in  the  region  of  interest,  the  peak  RMS  lift  coefficient  is  roughly  equal  for  the  flexible 
and  rigid  wings.  This  is  important,  as  it  signifies  that  the  lift  improvements  can  be  achieved  without  the  imposing 
high  lift  force  fluctuations.  The  level  of  buffet  decreases  significantly  with  sweep  angle,  with  peak  RMS  values  of 
the  60°  wing  reaching  just  a  third  of  that  observed  over  the  40°  wing. 

Figure  7  shows  the  variation  of  RMS  rolling  moment  coefficient  with  incidence  for  the  wings.  For  all  the  wings 
that  experienced  a  lift  enhancement,  a  significant  increase  in  RMS  rolling  moment  was  observed  in  the  lift 
enhancement  region  for  the  flexible  wing;  this  increase  was  not  observed  for  the  rigid  wings.  These  data  lend 
further  evidence  to  suggest  that  the  wings  vibrate  in  an  anti-symmetric  mode  in  the  lift-enhancement  region. 
Previously,  wing-tip  accelerometer  measurements  for  a  50°  wing  yielded  signals  from  each  side  of  the  wing  that 
were  180°  out-of-phase18.  In  the  current  case,  while  significant  rolling  moment  fluctuations  have  been  observed  in 
the  lift-enhancement  region,  the  fluctuations  in  lift  force  coefficient  were  not  significantly  larger  for  the  flexible 
wing  than  for  the  rigid. 

The  question  as  to  whether  or  not  an  anti-symmetric  vibration  of  the  wing  is  a  necessary  condition  for  the 
production  of  additional  lift  is  an  important  one,  and  one  that  is  as  yet  unanswered.  Flowever,  further  evidence  may 
be  presented  to  suggest  that  this  form  of  vibration  is  indeed  necessary.  Figure  8  shows  the  variation  of  wing-tip 
RMS  acceleration  with  incidence  for  a  half-wing  50°  flexible  model,  and  compares  this  with  the  full  wing  case  as 
presented  in  reference  18.  Although  direct  measurement  of  the  forces  over  the  half-wing  model  was  not  possible,  an 
increase  in  the  level  of  buffet  similar  to  that  shown  over  the  full-wing  model  would  indicate  the  existence  of  a  lift 
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enhancement  region.  However,  no  such  increase  in  buffet  was  detected  for  the  half-wing,  suggesting  that  the  wing 
did  not  experience  an  increase  in  lift. 

The  mode  shapes  and  natural  frequencies  for  the  wings  were  predicted  using  the  finite  element  analysis  software 
Ansys.  In  all  cases  the  mode  shapes  were  similar  to  those  documented  for  the  50°  wing  previously18,  but  the 
frequencies  at  which  each  mode  occurred  was  a  function  of  sweep  angle  as  demonstrated  in  Table  2. 


St 

A 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Mode  6 

40 

0.172  (S) 

0.172  (A) 

0.671  (A) 

0.683  (S) 

0.944  (A) 

0.955  (S) 

45 

0.201  (S) 

0.201  (A) 

0.738  (A) 

0.752  (S) 

1.064  (A) 

1.070  (S) 

50 

0.235  (S) 

0.235  (A) 

0.797  (A) 

0.811  (S) 

1.179  (S) 

1.220  (A) 

55 

0.276  (S) 

0.276  (A) 

0.862  (A) 

0.870  (S) 

1.069  (S) 

1.424  (A) 

60 

0.327  (S) 

0.327  (A) 

0.905  (S) 

0.932  (A) 

0.949  (S) 

1.679  (A) 

Table  2:  Finite  element  predictions  of  normalised  natural  frequencies  of  test  models 
(S  -  symmetric  mode,  A  -  anti-symmetric  mode). 


Spectral  analysis  of  the  rolling  moment  coefficient  was  undertaken,  and  dominant  peaks  were  observed  at 
frequencies  approximately  corresponding  to  the  2nd  anti-symmetric  modes  shown  in  Table  2.  These  data  are  plotted 
in  Figure  9,  whence  a  good  agreement  between  the  finite  element  prediction  and  the  experimental  results  can  be 
seen;  the  agreement  seems  to  improve  with  increasing  sweep  angle. 


B.  Flow  Topology 

An  initial  attempt  to  document  the  surface  flows  over  50°  rigid  and  flexible  wings  was  conducted  previously 
using  a  fluorescent  tuft  technique18.  It  was  shown  for  a  rigid  wing  at  a  =  21°  that  although  a  localized  region  of 
reattachment  was  observed  near  the  centerline,  the  dominating  features  of  the  flow  were  its  unsteady  nature  and  the 
existence  of  significant  regions  of  flow  reversal,  particularly  along  the  leading  edges.  For  the  flexible  wing, 
however,  a  much  more  energetic  flow  indicative  of  continued  reattachment  was  observed.  Although  flow  reversal 
was  still  a  feature  of  the  flow  towards  the  leading  edges,  the  flow  was  much  less  unsteady. 

The  disadvantage  of  the  tuft  visualization  technique  is  that  although  an  indication  of  flow  direction  may  be 
gleaned  from  the  video  footage  and  still  images,  no  indication  of  the  velocity  of  the  flow  at  any  point  is  provided. 
For  this  purpose  a  more  quantitative  method  of  measurement  was  required.  For  this  purpose,  PIV  measurements  of 
the  near-surface  flow  were  undertaken.  Figure  10  shows  the  magnitude  of  velocity  near  the  wing  surface  over  the 
rigid  50°  wing  at  a  =  21°.  Note  that  the  measurement  does  not  encompass  the  entire  wing  surface.  Figure  1 1  shows 
the  streamlines  calculated  from  the  velocity  vectors  of  the  same  dataset.  These  figures  confirm  that  the  flow  over 
the  rigid  wing  at  this  high  incidence  is  completely  stalled,  with  maximum  velocities  on  the  wing  surface  of  around 
u/Uoo  =  0.3  in  the  flow  reversal  region  near  the  leading  edges.  The  flow  towards  the  centerline  is  devoid  of  the  high 
axial  velocities  that  are  associated  with  shear  layer  reattachment. 

There  are  a  number  of  critical  points  that  may  appear  in  a  surface  flow  which,  upon  identification,  may  be  used 
to  interpret  the  three-dimensional  flow  topology.  A  critical  point  is  defined  as  a  point  in  the  flow  at  which  the  skin 
friction  and  the  velocity  vector  equals  zero20.  While  any  point  in  a  flow  can  have  only  a  single  trajectory,  the 
existence  of  a  singular  or  critical  point  in  a  flow  is  identified  by  the  convergence  of  skin-friction  lines  to  a  single 
point.  Examples  of  critical  points  are  nodes  (of  separation  and  attachment),  saddle  points,  and  foci;  all  of  which 
imply  a  different  characteristic  of  the  three-dimensional  flow.  Since  this  topic  has  been  covered  extensively  in  the 
literature,  and  a  detailed  discussion  falls  outside  the  remit  of  this  paper,  the  reader  is  directed  to  references  20  and  21 
for  further  reading. 

A  summary  of  the  critical  points  identified  over  the  rigid  50°  wing  is  shown  in  Figure  12a,  which  indicates  the 
existence  of  two  saddle  points  (Si  and  S2)  and  a  node  of  separation  (NO  up-stream  of  a  pair  of  foci  (F[  and  F2).  The 
total  number  of  nodes  must  exceed  the  total  number  of  saddles  on  an  isolated  body20  by  2,  but  clearly  in  the  current 
case,  since  the  foci  are  classified  as  nodes,  there  would  appear  to  be  a  node  missing  from  this  interpretation.  In  fact, 
since  we  are  considering  only  part  of  the  flow  this  is  acceptable,  provided  it  is  accepted  that  the  resulting 
interpretation  can  only  be  supposition  until  the  full  surface  flow  pattern  is  documented.  However,  the  measurements 
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have  elucidated  just  enough  of  the  surface  pattern  to  give  an  indication  of  the  separation  characteristics.  The  fact 
that  streamlines  are  directed  towards  a  pair  of  foci  in  the  manner  depicted  in  Figure  12  suggests  that  flow  leaves  the 
upper  surface  by  forming  a  pair  of  vortices,  as  shown  in  the  three-dimensional  interpretive  sketch  that  forms  Figure 
12a.  The  remainder  of  the  flow  is  directed  upstream,  and  the  source  of  this  is  likely  to  be  flow  ‘leaking’  over  the 
trailing  edge.  At  some  point,  due  to  the  pressure  gradient,  this  flow  forms  the  pair  of  vortices  that  curve  up  into  the 
free  stream,  while  the  remaining  surface  flow  continues  towards  the  apex,  where  it  will  separate  in  proximity  with 
the  separated  shear  layer. 

Figure  13  shows  the  magnitude  of  the  near  surface  time-averaged  velocity  over  the  flexible  wing  at  a  =  21°, 
while  Figure  14  shows  the  streamlines  calculated  from  the  velocity  vector  field.  Contrary  to  the  streamlines  over  the 
rigid  wing,  the  flexible  wing  demonstrates  a  symmetrical  time-averaged  flow.  Significant  regions  of  flow  reversal 
are  observed  over  the  wing,  with  high  negative  axial  velocities  peaking  at  u/U*,  =  -0.38.  Along  the  centerline,  high 
positive  axial  velocities  are  observed,  peaking  at  around  u/LL,  =  0.67.  A  node  exists  on  the  centerline  at  around  x/c 
=  0.2,  and  all  surface  streamlines  on  the  wing  issue  from  this  node,  indicating  that  it  is  a  node  of  attachment.  A  pair 
of  saddle  lines  exist  that  separate  the  flow  into  that  which,  having  attached  to  the  wing  surface,  continues  with  high 
velocity  directly  downstream,  and  that  which  becomes  entrained  into  the  region  of  tip  stall  and  flow  reversal  along 
the  leading  edges.  Figure  15a  shows  the  locations  of  the  identified  critical  points,  on  which  the  saddle  lines  are 
shown  as  (Si)  and  (S2).  Figure  15b  shows  the  three-dimensional  interpretation  of  the  time-averaged  flow.  This 
interpretation  has  the  shear  layer,  which  separates  from  the  leading  edges,  reattaching  inboard  of  the  leading  edges 
on  the  suction  surface,  and  curling  up  to  result  in  flow  reversal  in  this  region.  Flow  that  is  not  entrained  into  the  tip- 
stall  region,  passes  over  the  detached  shear  layer,  and  reattaches  aft  of  the  attachment  node  to  result  in  the  high  axial 
velocities  observed  along  the  centerline. 

Similar  streamline  patterns  have  also  been  observed  over  rigid  wings  of  38.7°  leading  edge  sweep22  at  a  =  10°, 
and  of  50°  sweep11  at  a  =  20°.  The  similarity  between  the  surface  flow  topologies  between  these  two  cases  is  very 
good,  despite  the  very  different  flow  conditions.  Flowever,  in  neither  case  was  an  attachment  node  recorded  along 
the  centerline,  indicating  that  the  flow  around  the  apex  does  not  separate  and  remains  attached  along  the  length  of 
the  wing.  This  is  likely  to  be  due  to  the  lower  incidences  of  these  wings  compared  with  the  flexible  wing  being 
studied  herein. 

The  data  therefore  show  that  the  flow  over  the  flexible  wing  is  much  more  coherent  than  that  over  the  rigid  wing 
at  the  same  incidence.  The  flexibility  clearly  results  in  fundamental  improvements  to  the  flow,  and  allows  the  shear 
layer  to  continue  to  reattach  to  much  higher  angles  of  attack  than  would  otherwise  be  expected. 

Off-surface  flows  were  also  documented  in  the  course  of  this  research.  Figure  16  shows  the  results  of  a  LDV 
measurements  of  the  chordwise  velocity  field  in  the  cross-flow  plane  at  x/c  =  0.7  for  the  (a)  rigid  and  (b)  flexible  50° 
wings  at  a  =  25°.  For  the  rigid  wing,  the  flow  is  much  more  symmetric  and  coherent  than  suggested  by  the  PIV 
measurements  presented  above  for  a  =  21°.  The  data  show  a  large  region  of  reversed  flow  extending  well  above  the 
surface  of  the  wing,  but  a  region  of  high  axial  velocity  remains  along  the  centre-line,  albeit  at  a  considerable 
distance  from  the  wing  surface.  Comparing  Figure  16  with  Figure  1 1  indicates  that  the  difference  of  2°  of  incidence 
between  the  two  measurements  was  significant.  At  a  =  25°,  the  LDV  measurements  indicate  that  the  flow  over  the 
rigid  wing  is  still  semi-attached,  while  at  27°  the  PIV  measurements  suggest  the  opposite  to  be  true. 

The  chordwise  velocity  field  over  the  flexible  wing  shares  many  similarities  with  that  of  the  rigid  wing.  Again,  a 
region  of  reversed  flow  extending  well  above  the  wing  surface  is  evident  along  with  a  region  of  high  axial  velocity 
above  the  surface  near  the  centre-line.  Flowever,  the  region  of  flow  reversal  is  reduced  in  extent  and  appears  flatter 
that  for  the  rigid  wing,  while  the  region  of  high  axial  velocity  above  the  centre-line  is  larger  and  extends  closer  to 
the  wing  surface. 

Figure  17  shows  the  axial  velocity  profiles  measured  above  the  centre-line  of  the  rigid  and  flexible  50°  wings  at 
(a)  x/c  =  0.4  and  (b)  x/c  =  0.7.  For  both  cases,  the  flexible  wings  exhibit  much  higher  axial  velocities  towards  the 
wing  surface  than  do  the  rigid  wings,  while  further  away  from  the  wing  the  velocity  profiles  merge. 


IV.  Conclusions 

Force  measurements  have  demonstrated  the  ability  of  a  flexible  non-slender  delta  wing  to  enhance  lift  and  delay 
stall  compared  with  a  rigid  wing  of  similar  geometry.  The  work  has  extended  the  results  of  a  recent  study  to  include 
data  for  a  range  of  wings  with  varying  sweep  angle.  Over  the  range  of  sweep  angles  tested  the  greatest  lift 
enhancement  was  observed  over  the  wing  with  smallest  leading  edge  sweep,  while  for  the  only  slender  wing 
considered  no  lift  enhancement  was  observed.  This  indicates  that  it  is  a  feature  of  the  fundamentally  different  flow 
that  occurs  over  low-sweep  wings  that  is  responsible  for  the  phenomenon. 
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Significantly,  the  variation  of  RMS  lift  force  with  incidence  for  the  rigid  and  flexible  cases  shows  good 
comparison,  while  the  RMS  rolling  moment  coefficients  for  the  flexible  wings  were  much  higher  than  the  rigid 
wings  in  the  lift  enhancement  region.  This  indicates  that  the  wings  are  vibrating  in  an  anti-symmetric  structural 
mode,  and  confirms  the  results  of  a  previous  investigation.  To  further  investigate  the  role  of  this  anti-symmetric 
vibration  in  the  lift  enhancement  phenomenon,  RMS  wing-tip  acceleration  measurements  for  a  full  and  half-wing 
50°  model  were  compared.  These  results  indicated  that  the  half-wing  model  did  not  experience  an  increase  in  lift, 
and  may  suggest  that  the  anti-symmetric  vibration  is  therefore  a  necessary  condition  for  the  lift  enhancement  to 
occur. 

PIV  and  LDV  measurements  have  demonstrated  the  striking  difference  between  the  flows  over  the  flexible  and 
rigid  wings  in  the  post-stall  region  at  a  =  25°  and  21°.  The  effect  of  flexibility  is  to  reduce  the  extent  to  the  region 
of  flow  reversal  that  occurs  over  the  wing  surface,  and  to  promote  the  reattachment  of  the  shear  layer  to  higher 
incidences. 
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Figure  1:  Schematic  of  the  working  section  and  high 
incidence  mechanism,  7'  x  5'  wind  tunnel. 
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Figure  3:  Principle  model  dimensions  of  flexible 
wing  with  leading  edge  sweep,  A  =  50°. 
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(a)  Locations  of  PIV  measurement  planes. 


Figure  2:  PIV  and  LDV  measurement  planes 
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(a)  A  =  40° 


(b)  A  =  50° 


(c)  A  =  60° 


Figure  5:  Variation  of  time-averaged  and  fluctuation  amplitude  of  wing-tip  deflections  with  incidence. 
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(a)  A  =  40° 
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Figure  7:  Variation  of  RMS  rolling  moment  coefficient  with  incidence. 
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Figure  8:  Variation  of  RMS  wing-tip  acceleration  for  50°  half-wing  model,  compared  with  that  of  full  wing 
model  from  reference  18. 


Figure  9:  Comparison  of  frequency  of  second  anti-symmetric  mode  predicted  by  finite  element  models  with 
dominant  frequencies  in  rolling  moment  spectra. 
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Figure  13:  Near  surface  time-averaged  velocity 
field  at  a  =  27°,  flexible  wing,  A  =  50°. 


x/c 


Figure  14:  Near  surface  time-averaged 
streamlines  at  a  =  27°,  flexible  wing,  A  =  50°. 


Figure  15:  (a)  2D  and  (b)  3D  topological  interpretations  of  the  near  surface  flow  over  the  flexible  50°  wing  at 
a  =  27°. 
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(a)  Rigid  Wing 


Figure  17:  Axial  velocity  profile  above  wing  centre-line  at  (left)  x/c  =  0.4  and  (right)  x/c  =  0.7  (A  =  50°  wing  at 
a  =  25°). 
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Physical  Mechanisms  of  Lift  Enhancement  for 
Flexible  Delta  Wings 

E.  Vardaki*,  I.  GursuP  and  G.  S.  Taylor* 

University  of  Bath,  Bath  BA2  7 AY,  United  Kingdom 


Passive  lift  enhancement  for  flexible  nonslender  delta  wings  has  been 
demonstrated  as  a  potential  method  for  the  control  of  vortex-dominated  wing  flows. 
Physical  mechanisms  of  lift  enhancement  and  the  effect  of  important  variables  are 
discussed.  Lift  enhancement  for  flexible  wings  with  low/moderate  sweep  is  a  very 
complex  phenomenon,  involving  self-induced  antisymmetric  vibrations  of  leading 
edges,  spanwise  camber  effect  due  to  the  large  time-averaged  deflection,  shear  layer 
instabilities,  reattachment,  increased  mean  vorticity  flux  and  circulation, 
re-formation  of  the  leading  edge  vortices,  possible  enhancement  of  streamwise 
pressure  gradient,  and  the  effects  of  frequency  and  edge  velocity. 


Nomenclature 


c 

f 

fc/Ux 

s 

t 

q 

CL 

E 

Rec 

S 

Ux 

a 

A 


p 

A 


chord  length 
frequency 

dimensionless  frequency 

semi-span 

wing  thickness 

free-stream  dynamic  pressure 

coefficient  of  lift 

elastic  modulus 

chord  Reynolds  number 

wing  reference  area 

free-stream  velocity 

angle  of  attack 

spanwise  bending  stiffness  normalized  by  free-stream  dynamic  pressure 
( =  12(l-v)pJJ Js3/Et3) 

Poisson’s  ratio 
density 

leading-edge  sweep  angle 


I.  Introduction 

PASSIVE  lift  enhancement  for  flexible  delta  wings  has  been  demonstrated  as  a  potential  method  for  the 
control  of  vortex-dominated  wing  flows1,2.  Force  measurements  over  a  range  of  non-slender  delta 
wings  (with  sweep  angles  A  =  40°  to  55°)  have  demonstrated  the  ability  of  a  flexible  wing  to  enhance  lift 
and  delay  stall  compared  with  a  rigid  wing  of  similar  geometry.  An  example  for  A  =  50°  is  shown  in 
Figure  1.  It  is  seen  that  lift  enhancement  was  achieved  in  the  post-stall  region.  Both  an  increase  in 
time-averaged  lift  coefficient  of  up  to  45%,  and  a  delay  in  stall  of  up  to  9  degrees  were  observed.  For  the 
only  slender  wing  considered,  of  60°  leading  edge  sweep,  no  lift  enhancement  was  observed,  indicating  that 
it  is  a  feature  of  the  fundamentally  different  flow  that  occurs  over  low  sweep  wings  that  is  responsible  for 
the  phenomenon.  The  main  objective  of  this  paper  is  to  review  the  physical  mechanisms  of  lift 
enhancement  and  to  discuss  the  most  important  variables  for  flexible  nonslender  delta  wings. 
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II.  Experimental  methods 

Wind  tunnel  experiments  were  conducted  in  the  high  speed  working  section  of  the  7x5  foot  closed 
circuit  facility  at  the  University  of  Bath.  Models  were  supported  using  the  high-incidence  rig  as  shown  in 
Figure  2,  and  various  leading  edge  sweep  angles  and  thickness  were  tested.  Each  model  was  a  planar  delta 
wing  with  a  pressure  surface  leading  edge  bevel  of  45°,  and  a  square  trailing  edge.  Experiments  were 
conducted  at  a  constant  free  stream  velocity  of  Uoo  =  31  m/s,  resulting  in  chord  Reynolds  numbers  in  the 
range  440,000  -  900,000.  Each  model  was  made  from  aluminium  alloy  1050A,  which  has  the  following 
material  properties:  p  =  2700  Kg/m3;  E  =  69  GPa;  and  v  =  0.3.  Maximum  blockage  for  the  wind  tunnel  was 
approximately  2.3%  for  the  60°  wing  at  a  =  50°.  Force  balance  measurements  were  undertaken  using  a 
six-component,  strain-gauged  internal  balance.  Forces  are  normalised  by  qS.  Surface  flow  visualizations 
were  achieved  using  lengths  of  thin  fluorescent  thread  taped  to  the  wing.  Quantitative  flow  measurements 
were  undertaken  using  Digital  Particle  Image  Velocimetry  (DPIV).  Illumination  of  the  plane  of  interest 
was  achieved  using  a  pair  of  pulsed  mini  Nd:YAG  lasers  with  a  maximum  energy  of  120  mJ  per  pulse. 
Images  were  captured  using  an  8-bit  digital  camera  with  a  resolution  of  4.2  million  pixels. 

In  addition  to  the  extensive  wind  tunnel  experiments,  complementary  experiments  were  conducted  in  a 
water  tunnel  facility  located  at  the  University  of  Bath.  The  tunnel  is  an  Eidetics  Model  1520  free  surface 
water  tunnel,  which  has  a  0.381  m  x  0.508  m  x  1.524  m  test  section  and  can  achieve  speeds  up  to  0.45  m/s 
through  a  closed  circuit  continuous  flow  system.  The  tunnel  has  four  viewing  windows,  three  surrounding 
the  working  section  and  one  downstream  allowing  axial  viewing.  The  height  of  the  test  section  above  the 
floor  is  adequate  to  allow  flow  visualisation  viewing  from  below  as  well  as  from  the  sides.  The  tunnel  also 
incorporates  dye  system  with  six  available  dye  tubes  to  enable  flow  visualization  with  different  colours. 

Experiments  were  conducted  using  an  aluminium  low  sweep  angle  delta  wing  model  with  50°  sweep 
angle.  The  model  had  a  45°  bevel  on  the  leading  and  trailing  edges  on  the  pressure  surface.  The  bevelling 
of  the  model  was  done  in  order  to  produce  a  sharp  leading  edge,  which  ensured  the  proper  formation  of  the 
leading  edge  vortices.  The  model  had  a  chord  length  of  c  =  89  mm  and  a  thickness  of  t  =  2  mm,  giving  a 
thickness  to  chord  ratio,  t/c  =  2.25%.  It  was  mounted  upside  down  on  the  tunnel  using  a  sting  projecting 
from  the  rear  of  the  model,  as  shown  in  Figure  2.  Additionally,  experiments  on  50°  sweep  angle  delta 
wings  with  various  flexibilities  were  carried  out.  Four  different  flexible  wings  were  tested  with  thickness 
of  0.1905  mm,  0.254  mm,  0.381  mm  and  0.508  mm  respectively.  The  former  two  wings  were  made  out  of 
polyester  and  the  latter  two  out  of  polypropylene. 

Visualisation  of  the  vortex  trajectories  was  achieved  using  food  colouring  dye,  diluted  1 :4  with  water.  A 
digital  video  camera,  with  a  capture  rate  of  25  frames  per  second  and  a  resolution  of  570,000  pixels,  was 
used  to  capture  images  from  the  dye  flow  visualisation,  and  was  interfaced  to  a  desktop  computer  via  the 
commercial  software  package  Pinnacle  Studio  DV,  enabling  real  time  viewing  of  the  wing  and  the  capture 
of  camera  images  and  video  recordings.  The  measurement  uncertainty  in  locating  the  vortex  breakdown 
position  was  approximately  2%  of  the  chord  length. 

For  PIV  measurements,  the  flow  was  seeded  with  commercially  available  hollow  glass  particles  of 
mean  diameter  of  4  pm  and  the  PIV  camera  was  placed  near  the  downstream  viewing  window,  as  shown  in 
Figure  3.  To  illuminate  the  desired  planes,  the  laser  system  was  placed  underneath  the  test  section  of  the 
water  tunnel,  as  shown  in  Figure  3.  A  combination  of  cylindrical  and  spherical  lenses  was  used  to  generate 
the  required  light  sheet.  Images  were  captured  using  an  8-bit  digital  camera  with  a  resolution  of 
2048  x  2048  pixels  and  a  maximum  capture  rate  of  7.5  frames  per  second,  producing  3.75  frames  in  cross¬ 
correlation.  The  commercial  software  package  Insight  v6.0  and  a  Flart  cross-correlation  algorithm  were 
used  to  analyse  the  images,  with  an  interrogation  window  size  of  32  by  32  pixels,  and  to  produce  velocity 
vectors  for  further  processing.  Sequences  of  30  instantaneous  frames  were  taken  for  each  case  and  the 
time-averaged  (or  phase-averaged  for  dynamic  roll  motion)  velocity  and  vorticity  fields  were  calculated. 
PIV  images  for  the  delta  wing  were  captured  at  80%  of  the  chord  length  and  at  a  zero  roll  angle  for  both  the 
dynamic  and  static  cases.  For  the  dynamic  case,  the  PIV  system  was  triggered  in  order  to  capture  the  flow 
field  at  zero  roll  angle  and  then  perform  phase  averaging.  The  Reynolds  number,  based  on  the  chord 
length,  was  Rec  =  26,000  for  all  the  experiments. 
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ITT.  Results 


A.  Antisymmetric  vibrations 

Vibrations  in  the  second  anti-symmetric  mode  (see  Figure  4a)  accompany  lift  enhancement.  Figure  4b 
shows  the  variation  of  root-mean-square  of  wing  tip  acceleration  as  a  function  of  incidence.  It  is  seen  that 
there  is  a  sudden  increase  of  tip  acceleration  in  the  region  corresponding  to  lift  enhancement.  These  large 
amplitude  vibrations  occur  as  a  result  of  flow/structure  interaction.  The  existence  of  self-induced 
antisymmetric  vibrations  in  this  region  is  similar  to  the  vortex  interactions  discussed  in  a  recent  review 
article3.  The  details  of  this  spectacular  flow/structure  interaction  remain  to  be  studied.  Also,  one  of  the 
remaining  unresolved  issues  is  the  role  of  the  antisymmetric  vibrations  and  whether  they  are  essential  for 
lift  enhancement.  Recent  evidence2  indicates  that  antisymmetric  interactions  are  essential  as  a  half-model 
delta  wing  did  not  exhibit  self-excited  vibrations  and  lift  enhancement. 

B.  Spanwise  camber 

Large  time-averaged  structural  deflections  also  exist,  which  produce  spanwise  camber.  Although  there 
are  significant  time-averaged  deflections  of  the  wing,  particularly  at  high  incidences,  these  are  unlikely  to 
contribute  to  the  lift  enhancement  phenomenon1.  Figure  5  shows  the  variation  of  lift  coefficient  with 
incidence  for  the  flexible  wing  and  a  rigid  curved  wing,  which  has  the  approximate  shape  of  the 
time-averaged  deflection  of  the  flexible  wing.  It  is  seen  that  the  curved  wing  does  not  exhibit  similar  lift 
enhancement. 

Note  that  no  streamwise  vortex  is  expected,  even  for  the  flexible  wing  in  the  lift-enhancement  region,  as 
the  leading  edge  vortices  already  break  down  at  the  apex  for  these  incidences.  However,  a  study  of  the 
effect  of  dihedral  applied  to  highly  swept  delta  wings  showed  an  overall  decrease  in  lift  coefficient4. 
Figures  3  and  5  suggest  that  this  effect  is  negligible  for  the  nonslender  wing  at  low  incidences  when  leading 
edge  vortices  exist.  Lift  force  prior  to  stall  is  unchanged,  while  it  is  only  in  the  post-stall  region  that  the  lift 
enhancement  is  observed. 

C.  Reynolds  number  effect 

In  general,  the  Reynolds  number  effect  is  more  pronounced  for  nonslender  wings5,  as  the  vortical 
structures  form  closer  to  the  wing  surface  and  strong  interactions  with  wing  boundary  layer  take  place. 
However,  even  for  a  nonslender  wing  with  sweep  angle  A  =  50°,  the  flow  approaches  an  asymptotic  state  at 
higher  Reynolds  numbers  (around  30,000),  with  further  increases  in  the  Reynolds  number  resulting  in  only 
small  variations  in  the  location  of  vortex  core  and  breakdown. 

A  comparison  of  time-averaged  velocity  in  a  plane  through  the  vortex  core  in  water  tunnel  experiments 
(Re  =  26,000)  and  in  wind  tunnel  experiments  (Re  =  620,000)  is  given  in  Figure  6.  It  is  seen  that  the 
time-averaged  breakdown  location  is  nearly  the  same,  although  there  are  some  differences  in  the  core 
structure  and  near  the  wing  centreline.  Furthermore,  since  the  lift  enhancement  is  observed  in  the  post-stall 
region  for  the  rigid  wing,  Reynolds  number  effect  is  expected  to  be  negligible. 

D.  Shear  layer  reattachment 

Surface  tuft  visualization  and  PIV  measurements  near  the  wing  surface  for  the  rigid  wing  show  that  the 
flow  is  fully  separated  in  the  post-stall  regime1’2.  However,  for  the  flexible  wing,  shear  layer  reattachment 
is  observed,  extending  the  region  of  partially  attached  flow  at  high  incidences.  The  main  mechanism  for  lift 
enhancement  seems  to  be  related  to  the  excitation  of  the  shear  layer  instabilities.  For  A  =  50°  wing,  the 
dominant  frequency  of  structural  vibration  was  around  St  =  0.75,  which  compares  well  with  the  dominant 
frequencies  of  the  shear  layer  instabilities6’7.  The  self-induced  vibrations  of  the  wing  energize  the  vortices 
shed  into  the  shear  layer  and  promote  reattachment  to  the  wing  surface,  resulting  in  delayed  stall  and 
increased  lift.  This  is  similar  to  the  response  of  the  flow  over  a  backward-facing  step  to  the  periodic 
excitation.  It  is  well  known  that,  for  both  laminar  and  turbulent  separation8,  excitation  enhances  the 
formation  of  vortical  structures  and  substantially  reduces  reattachment  length. 

Flow  visualizations  for  a  rigid  and  various  flexible  wings  are  shown  in  Figure  7  for  a  =  20°  and 
A  =  50°.  Here,  the  parameter  As  is  the  spanwise  bending  stiffness  normalized  by  free-stream  dynamic 
pressure,  i.e.,  As  =  1 2(  1  -v2)pxU,.2s3/Et3,  where  E  is  elastic  modulus,  v  Poisson’s  ratio,  s  semi-span,  and 
t  thickness  of  the  wing.  The  dye-free  region  near  the  wing  centreline  is  much  broader  for  the  flexible 
wings,  indicating  earlier  reattachment  of  the  separated  flow.  For  the  most  flexible  wing  (/Ls  =  21.28),  it  is 
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interesting  that  re-formation  of  the  leading  edge  vortex  and  rapid  breakdown  are  visible  in  a  region  very 
close  to  the  wing  apex.  This  will  be  discussed  in  detail  later  on  in  the  paper. 

Figure  8  shows  the  vectors  and  contours  of  magnitude  of  time-averaged  velocity  in  a  crossflow  plane  at 
x/c  =  0.80  for  rigid  and  flexible  wings  for  an  even  larger  incidence  of  a  =  25°.  For  the  rigid  wing,  there  is 
no  reattachment,  however,  for  the  flexible  wings,  there  is  larger  velocity  near  the  wing  centerline  and 
reattachment  in  the  time-averaged  sense. 

E.  Mean  vorticity  flux 

For  the  rigid  and  flexible  wings,  the  vorticity  distributions  in  a  crossflow  plane  at  x/c  =  0.80  based  on 
PIV  measurements  are  shown  in  Figure  9  for  a  =  20°.  There  are  substantial  effects  on  the  vortical  flow 
with  increasing  wing  flexibility,  as  the  oscillating  leading  edge  is  an  unsteady  source  of  vorticity.  In  fact, 
the  time-averaged  vorticity  flux  is  expected  to  increase  with  the  vibrations  of  the  leading-edge,  as  this  term 

is  proportional  to  U ^  ,  where  Us  denotes  the  velocity  outside  the  boundary  layer  at  the  separation  point. 
This  is  confirmed  by  the  larger  values  of  vorticity  for  the  flexible  wings  in  Figure  9.  The  variation  of 
circulation  with  the  spanwise  bending  stiffness  normalized  by  free-stream  dynamic  pressure,  is  shown  in 
Figure  10.  It  is  seen  that  the  circulation  of  the  vortical  flow  in  a  crossflow  plane  increases  with  flexibility 
in  comparison  to  the  rigid  wing. 

F.  Vortex  re-formation 

In  order  to  simulate  the  effect  of  antisymmetric  vibrations,  experiments  with  a  rigid  wing  undergoing 
small  amplitude  rolling  motion  were  also  conducted.  This  approach  has  the  advantage  of  independent 
control  of  the  frequency  parameter,  which  cannot  be  varied  for  a  flexible  wing  without  affecting  the 
amplitude  of  vibrations.  Figure  11  shows  flow  visualization  for  the  stationary  and  rolling  delta  wing  (with 
an  amplitude  of  5°)  for  an  incidence  of  a  =  25°.  It  is  seen  that  the  totally  separated  flow  for  the  stationary 
wing  becomes  very  organized  for  the  rolling  wing  with  increasing  frequency.  The  most  interesting 
observation  is  the  re-formation  of  the  leading  edge  vortices  at  high  frequencies.  Although  the  leading  edge 
vortices  become  stronger  due  to  the  leading  edge  motion,  vortex  breakdown  is  delayed  for  the  rolling  wing 
compared  to  the  stationary  wing  for  which  breakdown  is  at  the  apex.  This  appears  to  be  in  contrast  to  the 
well-known  studies  of  vortex  breakdown,  which  indicate  that  increased  strength  of  vortices  should  cause 
premature,  rather  than  delayed,  breakdown.  This  result  suggests  that  streamwise  pressure  gradient  might  be 
modified  favorably  due  to  the  wing  motion. 

G.  Optimum  frequency 

Figure  12  shows  the  variation  of  mean  breakdown  location  as  a  function  of  dimensionless  frequency  for 
different  angles  of  attack  in  the  range  a  =  20°  to  30°.  It  is  seen  that,  for  each  incidence  in  the  post-stall 
region,  the  breakdown  location  is  zero  for  the  stationary  wing  (fc/Uoo  =  0).  Maximum  delay  of  the  vortex 
breakdown  location  is  achieved  in  the  range  of  fc/Uoo  —  1  to  2.  Again,  this  range  compares  well  with  the 
dominant  frequencies  of  the  shear  layer  instabilities6'7  for  a  nonslender  wing  of  A  =  50°. 

Figure  13  compares  the  variation  of  mean  breakdown  location  as  a  function  of  dimensionless  frequency 
for  two  amplitudes  of  rolling  motion  for  a  =  25°.  Even  for  the  small  amplitude  of  1°,  it  is  possible  to  have 
re-established  leading  edge  vortices.  The  optimum  dimensionless  frequency  is  also  in  the  same  range. 

H.  Leading  edge  velocity 

In  the  case  of  the  rigid  wing  undergoing  small  amplitude  rolling  motion,  the  dimensionless  frequency 
fc/Uoo  is  not  only  a  ratio  of  time  scales  of  convective  time  and  rolling  motion,  but  also  the  ratio  of  the 
leading  edge  velocity  and  freestream  velocity  for  a  given  wing.  Figures  12  and  13  suggest  that  increasing 
velocity  of  leading  edge  (with  increasing  frequency  parameter)  does  not  necessarily  delay  breakdown 
location.  Flowever,  it  is  expected  that  the  time-averaged  vorticity  flux  will  increase  with  increasing 
velocity  of  the  leading  edge.  This  is  confirmed  by  the  vorticity  distributions  in  a  cross-flow  plane  at 
x/c  =  0.80  as  shown  in  Figure  14  for  a  =  25°.  Note  that  the  results  shown  in  this  figure  are  phase-averaged 
over  30  cycles,  and  correspond  to  zero  roll  angle  as  the  roll  angle  changes  the  sign  from  negative  to 
positive.  This  is  the  reason  for  slight  asymmetry,  which  is  the  result  of  the  well-known  hysteresis  effect.  It 
is  seen  that  larger  values  of  vorticity  exist  with  increasing  frequency.  The  variation  of  circulation  with  the 
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dimensionless  frequency  is  shown  in  Figure  15.  It  is  seen  that  the  circulation  of  the  vortical  flow  in  a 
crossflow  plane  increases  with  frequency  in  comparison  to  the  stationary  wing. 


IV.  Conclusions 

Passive  lift  enhancement  for  flexible  nonslender  delta  wings  has  been  demonstrated  as  a  potential 
method  for  the  control  of  vortex  dominated  wing  flows.  Physical  mechanisms  of  lift  enhancement  and  the 
effect  of  important  variables  have  been  discussed.  Lift  enhancement  for  flexible  wings  with  low/moderate 
sweep  is  a  very  complex  phenomenon.  The  generation  of  self-induced  antisymmetric  vibrations  of  leading 
edges  is  not  completely  understood,  however  it  is  certain  that  these  vibrations  produce  the  observed  lift 
enhancement.  The  spanwise  camber  effect  due  to  the  large  time-averaged  deflection  does  not  contribute  at 
all.  The  main  mechanism  of  lift  enhancement  is  the  excitation  of  shear  layer  instabilities  and  promotion  of 
reattachment  of  the  separated  flow.  There  are  substantial  effects  on  the  vortical  flow  with  increasing  wing 
flexibility,  as  the  oscillating  leading  edge  is  an  unsteady  source  of  vorticity.  The  time-averaged  vorticity 
flux  increases  due  to  the  oscillating  leading  edge,  which  leads  to  increased  circulation.  Despite  the 
increasing  strength  of  the  vortical  flow,  re-formation  of  the  leading  edge  vortices  with  axial  flow  is 
observed,  which  suggests  that  the  streamwise  pressure  gradient  might  be  modified  favorably  due  to  the 
wing  motion.  The  mean  breakdown  location  becomes  a  maximum  at  an  optimum  frequency. 
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Plate 


Drive  Motor 


Figure  1:  Variation  of  lift  coefficient  with 
incidence  in  wind  tunnel  experiments  for  a  50°- 
sweep  angle  delta  wing. 


Figure  2:  Schematic  of  the  working  section  and 
high  incidence  rig,  T  x  5’  wind  tunnel,  University 
of  Bath. 


Axial 


Figure  3:  Overview  of  experimental 
setup  in  the  water  tunnel. 


(a) 


Figure  4:  (a)  Dominant  anti-symmetric 
mode  of  vibration  in  wind  tunnel 
experiments;  (b)  variation  of  root-mean- 
square  of  wing  tip  acceleration  as  a  function 
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Figure  8:  Vectors  and  contours  of  magnitude  of  time-averaged  velocity  in  a  cross-flow  plane  at 
x/c=0.8  for  rigid  and  flexible  wings  in  water  tunnel  experiments,  a  =  25°. 
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Figure  9:  Magnitude  of  vorticity  in  a  cross-flow  plane  at  x/c=0.8  for  rigid  and  flexible  wings  in  water 
tunnel  experiments,  a  =  20°. 
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Figure  10:  Variation  of  normalized  circulation  with  wing  flexibility. 
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Figure  11:  Flow  visualization  for  a  stationary  and  small-amplitude  (4>=50)  rolling  wing  in  water 
tunnel  experiments. 
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a  =  25  deg,  <)>  =  5  deg 


Figure  13:  Variation  of  mean  breakdown 
location  as  a  function  of  dimensionless 
frequency  for  two  values  of  rolling  motion  in 
water  tunnel  experiments,  a=25°. 


counter-clockwise  vortex 


Figure  15:  Variation  of  normalized  circulation  of 
vortical  flow  in  a  cross-flow  plane  at  x/c  =  0.8  as  a 
function  of  dimensionless  frequency  in  water 
tunnel  experiments,  a=25°. 
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Figure  14:  Magnitude  of  vorticity  in  a  cross-flow 
plane  at  x/c=0.8  for  stationary  and  rolling  wings  in 
water  tunnel  experiments,  a  =  25°. 
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